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Abstract 

Motivated by recent discussions of the string-theory landscape, we propose 
field-theoretic realizations of models with large numbers of vacua. These mod- 
els contain multiple ?7(1) gauge groups, and can be interpreted as deconstructed 
versions of higher-dimensional gauge theory models with fluxes in the compact 
space. We find that the vacuum structure of these models is very rich, defined 
by parameter-space regions with different classes of stable vacua separated by 
boundaries. This allows us to explicitly calculate physical quantities such as 
the supersymmetry-breaking scale, the presence or absence of i?-symmetries, 
and probabilities of stable versus unstable vacua. Furthermore, we find that 
this landscape picture evolves with energy, allowing vacua to undergo phase 
transitions as they cross the boundaries between different regions in the land- 
scape. We also demonstrate that supergravity effects are crucial in order to 
stabilize most of these vacua, and in order to allow the possibility of cancelling 
the cosmological constant. 
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1 Introduction 



Recent developments in the study of string-theory compactifications suggest that 
there exist huge numbers of string vacua, with different cosmological constants and 
different low-energy phenomenological properties [1]. The resulting picture, dubbed 
the "landscape" [2], has stimulated a statistical analysis of the number of string 
vacua [3], the supersymmetry- breaking scale, and other phenomenological features [4]. 
Anthropic principles have even been advanced to resolve difficult issues such as the 
cosmological constant problem [5]. 

One of the problems facing this landscape picture of string theory is that of calcu- 
lating physical parameters. This is, to a large extent, due to the limited technology 
for performing string calculations in the relevant flux vacua [6]. It is therefore use- 
ful to present field-theoretic counterparts of such constructions, i.e., field-theoretic 
models which naturally give rise to very large numbers of vacua, and to be able to 
quantitatively determine statistical distributions of relevant physical quantities such 
as the cosmological constant, the scale of supersymmetry breaking, the Higgs mass, 
gauge and Yukawa couplings, and the like. 

The purpose of this paper is to present such field-theoretic examples based on 
multiple Abelian gauge groups and multiple charged scalar fields. As wc shall see, 
such models naturally lead to large numbers of vacua, and allow us to quantitatively 
address many of the pressing questions that such pictures raise. 

In Sect. 2, as a brief warm-up, we discuss the standard Fayet-Iliopoulos model, 
now reconsidered from our perspective. In Sect. 3, we then describe, in some detail, 
the construction of a relatively simple toy model which nevertheless gives rise to 
many of the landscape features that we wish to study. Specifically, we consider a su- 
persymmetric model with two C/(l)'s with Fayet-Iliopoulos terms, and three charged 
chiral superfields. We show that this model contains seven distinct classes of vacua 
which populate the parameter space in different pie-slice regions, giving rise to vastly 
different physical properties such as supersymmetry-breaking scale and numbers of 
unbroken gauge and global symmetries. We also show that renormalization-group 
evolution can bring the theory from one vacuum to another, while the boundaries 
between different regions correspond to parameter-space surfaces on which a single 
scalar particle becomes massless. 

In Sect. 4, we consider a generalization of this toy model to more complicated 
supersymmetric models with large numbers of U{1) gauge groups. Such models eas- 
ily give rise to huge numbers of distinct classes of vacua. The simplest example is 
provided by n copies of the Fayet-Iliopoulos model, thereby giving rise to 3" distinct 
classes of vacua. We concentrate most of our efforts, however, on a four-dimensional 
model which we believe is the "deconstruction" [7] of a six-dimensional model with 
internal magnetic fluxes. Our interest in this respect is to provide a purely four- 
dimensional description of these higher-dimensional flux models which are among the 
simplest models that generate chirality in four dimensions. This flux-compactiflcation 
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framework is also probably the simplest and most successful framework for the con- 
struction of semi- realistic string models with intersecting branes [8] . We also comment 
on the role of mixed anomalies and their cancellation in the vacua we consider, and 
show that almost all of these vacua are unstable, in analogy with the Nielsen-Olesen 
instabilities of internal magnetic-flux compactifications. 

In Sect. 5, we demonstrate that supergravity effects can lead to large soft mass 
terms which stabilize most of these unstable vacua in certain regions of the parameter 
space. We also show that qualitative arguments based on cosmological constant 
cancellation generically require a large gravitino mass which, in a gravity-mediated 
supersymmetry breaking scenario, is consistent with the large soft terms that stabilize 
our unstable vacua. Finally, we conclude with comments about possible applications 
and generalizations of these constructions, and prospects for future research. 

2 Starting small: The Fayet-Iliopoulos model 

The simplest example of a model exhibiting a non-trivial, multiple- vacuum struc- 
ture that we have in mind is the Fayet-lliopoulos (FI) model of supersymmetry break- 
ing [9]. Recall that this model has a single C/(l) gauge group with gauge couphng g 
and FI coefficient ^, two charged chiral superfields of charges ±1, and a super- 
potential W = m$'^+)$(~^. It turns out that the scalar potential for this model has 
three classes of extrema: 

• Extrema with v^^^ = {(f)^^^) = 0. Such extremum solutions exist for all ^, but 
represent stable vacua only for ± g^^ > 0. 

• Extrema with ^ 0, v^~'> — 0. These solutions exist if -|- g'^^ < 0, and 
are stable whenever they exist. 

• Extrema with v'^"^ ^ 0, v*^"*"^ = 0. These solutions exist if vr? — g^^ < 0, and 
are stable whenever they exist. 

Note that in the first class of solutions, both the gauge symmetry and /^-symmetry 
of the model are unbroken, whereas in the second two classes of solutions they are 
both broken.* This vacuum structure is illustrated in Fig. 1, which indicates the 

appropriate vacuum as a function of 

Of course, the FI model is actually defined by more than the single parameter 
^; we can just as easily consider m as an additional defining parameter. We would 

*Z?-term breaking always preserves i?-symmetries, independently of their detailed definition. By 
contrast, F-type breaking only breaks -R-symmetries if the auxiliary field F is charged under the R- 
symmctry. In all of the examples to be considered in this paper, wc find that _F-tcrms arc charged 
under the i?-symmetry if we define i?-syminetry by requiring the superpotential to be invariant 
and assume that all matter fields have the same i?-charge. Thus, F-breaking is equivalent to R- 
symmctry breaking. Indeed, for the FI model, if the superpotential charge is equal to 2, then the 
matter superfield i?-charges are 1 and therefore the F-terms have charges —1. 
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Figure 1: The vacuum structure of the Fayet-Ihopoulos (FI) model as a function of the 
FI coefficient ^. The different classes of vacua outlined in the text (and here denoted 
{0}, {+}, and {— } respectively) occupy non-overlapping regions along the ^-axis. In 
the regions corresponding to {+} and {— } vacua, we have indicated the existence of an 
unstable {0} extremum in parenthesis. We have also shaded the regions corresponding 
to vacua with at least one non-zero vacuum expectation value, corresponding to vacua 
exhibiting i2-symmetry breaking. 




Figure 2: The vacuum structure of the FI model, now sketched as a function of (,^,m). 
Each region now has the shape of a pie-slice. The slopes of the diagonal lines are iy^. 

then sketch the vacuum structure of this model as in Fig. 2, which shows the vacuum 
structure as a function of both ^ and m. The slopes of the diagonal hnes are ±5^^. 

Finally, we may recognize that g also represents an additional degree of freedom 
defining the model. In this case, it is relatively easy to visualize the resulting vacuum 
structure as a function of location in a three-dimensional space. 

Suppose we were to ask a simple question: given the vacuum structure of this 
model, what is the net probability of randomly finding a vacuum that preserves R- 
symmetry? In this simple example, the answer clearly depends on the variables over 
which we average. If we consider ^ to be the only variable as in Fig. 1 and perform 
our averaging with measure d^^ we would claim that we have a virtual impossibility 
of realizing such a vacuum for any given (fixed) value of {mj gY . This is because the 
{0} region in Fig. 1 ultimately occupies only a region of finite size (set of measure 
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P = - I d9jl--9,] = -. (2.3) 



zero) within the infinitely long ^-axis.^ On the other hand, if we consider both ^ and 
m to be our relevant variables as in Fig. 2 and calculate our average with measure 
d^dm^ — 2'md^dm, we would find that this probability is given by 

2 

Prob (R symm) = 1 tan~^ . (2-1) 

TT 

Indeed, this is nothing but the ratio of the angles of the pie-shces in Fig. 2. Of 
course, this result would change if we were to take our measure to be, for example, 
d^dm rather than d^dm^. Finally, if we wish to consider g as a, third defining degree 
of freedom, we need to integrate the result in Eq. (2.1) over possible values of g. 
Clearly, if we use the measure dg, we would obtain a vanishing probability because 
of the infinite relative weighting given to vacua with large g. However, motivated by 
the pie-slice diagram in Fig. 2, we can define an angular variable 9g = tan~^ g^ and 
take our measure to be 

This naturally suppresses the contributions from laige-g vacua, and yields the inte- 
grated result 

fV2 / 2 \ 1 
= 2 

Of course, in practice we should perhaps only be integrating up to the limits of our 
perturbative treatment, restricting to g ^ y/An, but this would introduce a possible 
strong dependence on the cutoff employed. Likewise, the measure in Eq. (2.2) is only 
one of many possible measures we could use. In any case, this discussion dramat- 
ically illustrates the importance of clearly defining both our relevant variables, our 
integration measure, and all associated cutoffs when averaging over multiple vacua. 

It may seem strange to consider a coupling constant as a coordinate parametrizing 
a space of vacua. However, we recall that in string theory, the coupling constant is 
related to the vacuum expectation value (vev) of the dilaton, and such a vacuum 
expectation value is traditionally taken as one coordinate of an ensuing moduli space 
of vacua. 

The issue of determining the appropriate landscape coordinates is clearly im- 
portant, but unfortunately model-dependent. In the popular compactifications with 
fiuxes [6], the gauge couplings on D3-branes depend on the dilaton which, in turn, 

depends on the fluxes. However, since we ultimately have to restrict our attention 
to gauge couplings at high scales which are compatible with the experimentally mea- 
sured couplings at low scales, these couplings should be ~ 0{1). A pragmatic attitude 

'''In practice, of course, we must have ^ < Mp where Mp is the Planck scale, so our parameter 
space is not truly infinite in extent. Also, the FI coefficients ^ are typically quantized. Throughout 
this paper, however, we shall consider the continuous limit and imagine that our parameter spaces 
are infinite in extent. 
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would then be to focus directly on that subset of fluxes leading to realistic gauge cou- 
plings. By contrast, the FI terms depend on moduli fields and charge assignments; 
they are certainly landscape-dependent and are therefore natural landscape coordi- 
nates. The super symmetric mass in the FI model can also be generated by suitable 
fluxes, in analogy with the /x-term of MSSM [10]. Likewise, the supergravity-induced 
soft mass terms that we shall introduce in subsequent sections are also flux-dependent. 
Thus, we see that essentially all mass-like parameters are natural landscape coordi- 
nates. 

An important constraint is that the FI terms and supersymmetric and soft mass 
terms are bounded in any string compactifications. Indeed, FI terms are proportional 
to vev's of various moduh fields and to the sum of U{1) charges, and are always smaller 
than the string and the Planck scale. Likewise, the various fluxes must generically 
satisfy Ramond-Ramond tadpole constraints; in most cases, these constraints place 
upper limits on their possible values. Thus, in Figs. 1 and 2 (as well as all subsequent 
parameter-space diagrams in this paper) , we should implicitly impose a cutoff on the 
FI terms and various mass parameters. This cutoff is model-dependent, but should 
roughly be of size ~ 0{Mp) where Mp is the Planck mass. 

An interesting result of our analysis of the FI model in Fig. 2 is that the vacuum 
structure is energy- dependent. Indeed, if Z„i{E) denotes the wavefunction renormal- 
ization factor for the matter fields (which we take to be equal for the two superfields) 
and Z{E) denotes the wavefunction renormalization factor for the gauge superfield, 
then the boundaries between the {-|-} and {0} regions and the {0} and {— } regions 
will scale according to the energy-dependent expressions ±g^Zjn(E)/Z(E). These 
boundaries therefore rotate under the renormalization group. This opens up the pos- 
sibility of having a vacuum which passes from one region to another under RG flow, a 
phenomenon which can generate phase transitions with potentially interesting cosmo- 
logical implications. We will study this phenomenon in more detail in Sect. 3, within 
the framework of a more interesting (and also slightly more complicated) example. 

3 Constructing a toy model 

It turns out that the FI model is not sufficiently complex for our purposes. In this 
section, we shall therefore turn to the next-simplest model which, as we shall see, will 
exhibit almost all of the relevant features which we shall encounter when we proceed 
to consider more complicated situations in subsequent sections. In particular, we shall 
see that this toy model gives rise to a non-trivial "landscape" consisting of multiple 
stable (or metastable) vacua with different low-energy phenomenologies, unstable 
extrema, phase transitions, etc. Moreover, even though this toy model is relatively 
simple, we expect that the resulting landscape is literally a component of the full 
string-theory landscape in cases of string models with multiple U{1) gauge factors. 

Our model consists of two U{1) gauge symmetries, denoted f/(l)i and U{1)2, and 
three J\f = 1 chiral superfields, $1=1,2,3- The charges of these chiral superfields under 
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$1 


-1 


$2 


1 -1 


$3 


1 



Table 1: C/(l) charge assignments for chiral superfields in our toy model. 



the C/(l) gauge symmetries are chosen as in Table 1. In a string-theory context, such 
U{1) gauge factors can be imagined as arising from different D-branes, and the 
fields can arise as strings stretched between these branes. We shall also assume that 
the J\f = 1 supersymmetry is broken by Fayct-Iliopoulos D-term coefficients and 
^2, and by a renormalizable Wilson-line supcrpotential of the form 



(3.1) 



Our model is thus defined by three parameters, {'Ci)'C2)^}) and our goal will be to 
study the vacuum structure of this model as a function of these parameters. Of 
course, the resulting physics is unchanged if A — > —A. We shall therefore restrict 
ourselves to A > for simplicity. 

It is straightforward to analyze the vacuum states of this theory. As usual, the 
scalar potential is given by 



a 

where the D- and F-terms are given by 



Em 



dW 

d(pi 



(3.2) 



(3.3) 



Here a = 1,2 is the index of the gauge group U{1) factor, Qa is the gauge coupling 
corresponding to the U{l)a factor, and i = 1, 2, 3 is the index of the chiral supcrficld. 
Thus, q^""^ denotes the U{l)a charge of $i. In most of our considerations the gauge 
couplings Qa will not be important, so we will henceforth consider gi — g2 — 1 
for simplicity. We will reinstall gauge couplings whenever relevant for our analysis. 
Our task is to determine the extrema of V by seeking solutions to the simultaneous 
equations 

dV dV 

^ = ^ = ' (3.4) 
d(f)i d(f)* 

and then to determine whether these extrema represent stable (or metastable) vacua 
by calculating the eigenvalues of the corresponding mass matrix 



(3.5) 
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evaluated at the extrema. Note that in general, there will be a zero eigenvalue 
for each spontaneously broken C/(l); these eigenvalues correspond to the resulting 
Nambu-Goldstone bosons. The extrema defined by Eq. (3.4) represent stable (or 
metastable) vacua only if each of the remaining eigenvalues is positive. 

Just as with the FI model in Sect. 2, it will prove convenient to group the extrema 
of V into classes depending on which combinations of chiral superfields receive non- 
zero vacuum expectation values. This classification will help in determining such 
features as whether the sources of SUSY-breaking arc primarily D-terms or F-tcrms, 
and whether they are likely to lead to /^-symmetry breaking when incorporated into 
a supergravity framework. Since there are three chiral superfields in this example, 
there are correspondingly 8 = 2^ classes of extrema of V. Denoting = Vj, we 
shall define our classes of extrema according to their values of Vi, using the notation 
{ijk...} to indicate the class of vacua in which Vi, Vj, Vk, ... are all non-zero (with 
{0} denoting the vacua in which all Vi vanish). Note that in this toy model, we can 
choose all Vi to be real without loss of generality. 

3.1 A = l 

Let us first consider the vacuum structure of this model when A = 1. Our results 
are then as follows. We find that {0} extrema exist for all values of (^^1,^), but are 
always unstable. This result, of course, is intuitively obvious, since {0} vacua are 
essentially at the unstable center of the SUSY-breaking "Mexican hat". By contrast, 
{1} extrema exist only for > 0. However, these extrema are stable only within 
the subregion defined by |^2| < Ci- Likewise, {2} extrema exist for ,^2 > Ci ^^'^ ^-^e 
stable in the region {1 < 0, ^2 > 0, while {3} extrema exist only for ^2 < and 
are stable within the subregion |^i| < |^2|- Similarly, {12} extrema exist only for 
< ^1 < ^2 and arc stable everywhere within this region, while {23} extrema exist 
only for ^1 < ^2 < and are also stable everywhere within this region. Finally, for 
A = 1, extrema in Classes {13} and {123} do not exist for any values of (^1,^2)- 

It is clear from these results that for A = 1, there are no regions in {^1,^2) 
parameter space in which there are multiple stable vacua. Instead, the different 
classes of vacua occupy non-overlapping regions which completely fill out this space. 
This situation is sketched in Fig. 3. 

It is important to interpret the diagram in Fig. 3 correctly. Each point in the 
(^1,^2) plane represents a unique theory, i.e., a unique model parametrized by its 
defining parameters ^1 and ^2- Thus, this plane is nothing but a field-theoretic 
"landscape" for this two-C/(l) theory with three chiral superfields. For each point 
in the landscape, there is a unique vacuum whose energy can be calculated from 
Eq. (3.2). For example, for models within the region specified by ^2 > —^1 and 
^2 < {i.e., the lower half of the {1} region in Fig. 3), the energy of the stable {1} 
extremum is given by V{iy — ^|/2. This is to be contrasted with the unstable {0} 
and {3} extrema which also exist for models in this region, with energies given by 
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{0}{1}{2} 




^1 



{0}{3} 



Figure 3: The "landscape" of our toy model with two U{1) gauge fields and three chiral 
superfields, sketched for A = 1. The different classes of vacua (labelled according to which 
chiral fields receive non-zero vacuum expectation values) occupy the non-overlapping regions 
indicated above in bold type, while unstable extrema in each region are indicated in smaller 
type. We have shaded the regions corresponding to vacua with two non-zero vacuum 
expectation values (Classes {12} and {23}); these are regions in which i?-symmetry is 
potentially broken. 



V{0} = + ^|)/2 and V{3} = respectively. Note that indeed V{i} < V{0},{3}. 

At this stage, we can already reach an interesting conclusion: for the A = 1 
landscape, the average ratio of stable extrema relative to all extrema is only 

, , n n , n 1/11111111 

Probabihty that extremum is stable = - -H 1 1 1 1 1 1 — 

^ 8 V4 2332422 

19 

= — fti 0.396 . (3.6) 
48 ^ ' 

Note that in calculating this probability, we are averaging over the landscape with a 
uniform measure (averaging the stability probabilities in each of the eight octants of 
the plane), imagining that all values of our landscape coordinates (^1,^2) are equally 
likely. This result is significantly larger than the value (1/2)^ which we might have 
naively expected from the fact that our mass matrix in Eq. (3.5) has six eigenvalues, 
any one of which randomly could have been negative. Even if we allow for potentially 
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as many as two zero eigenvalues corresponding to Nambu-Goldstone bosons for the 
two U{1) gauge factors in this toy model, we still have a much higher probability of 
finding stable extrema than expected. 

In a similar vein, we can also see from Fig. 3 that a randomly chosen model in 
the landscape has exactly a 25% chance of breaking ii-symmetry when A = 1. Those 
regions populated by i?-symmetry breaking vacua are shaded in Fig. 3. 

In this landscape, there do not exist any regions with overlapping stable vacua. In 
other words, there are no values of (^i, ^2) giving rise to multiple stable vacua. Thus, 
each vacuum is absolutely stable with respect to dynamics within that model. How- 
ever, these different vacua fill out a landscape in theory space, as sketched in Fig. 3, 
and these vacua ultimately have different energies which depend on the landscape 
coordinates (^1,^2)- For A = 1, the different classes of extrema have energies: 



Moreover, note from Eq. (3.7) that Class {1} vacua with {2 = 0, Class {2} 
vacua with ^2 = and Class {3} vacua with ,^1 = are all SUSY-preserving, 

with V = 0. Even though these SUSY-preserving vacua only form a one-dimensional 
space of vacua within the larger two-dimensional landscapes we have plotted (thereby 
essentially forming a set of measure zero within the full landscape) , these vacua rep- 
resent the lowest-energy vacua within their respective regions. However, unless tran- 
sitions between difi'crent vacua in the landscape are somehow possible, the existence 
of such super symmetric vacua at one landscape location does not destabilize the 
non-supersymmetric vacua that may exist at other landscape locations. 

We see, then, that this "landscape" is indeed a very realistic, and ultimately 
completely calculable, model of the string-theoretic landscape. Indeed, we have a 
space of models, with each parametrized by the model-defining parameters (^1,^2), 
and each with a different vacuum and vacuum energy. 

3.2 A < 1 

It is important to consider how this landscape picture is modified for A 7^ 1. 
Considering first the case with A < 1, we find that Class {0} extrema continue to 
exist for all (^i,S) ^-^d continue to be unstable everywhere. Class {1} extrema, by 
contrast, continue to exist for all .^1 > 0, but are now stable only within the smaller 
region |^2| < ^^Ci- Similarly, Class {2} extrema continue to exist for all ^2 > 



^{0} 

V{12} 
V{23} 



|(e? + ^2') = ^{l} + ^{3} 

2 '^2 

i (6 + 6)' 

lf2 

+ 266+^2) = V|2}-V|3} 



(3.7) 
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are stable only within the region defined by 



(1 - A^) 6 < - (A^ + 1) 6 , (A^ + 1) 6 > (A^ - 1) 6 , (3.8) 

while Class {3} extrema continue to exist for all S < but are stable only within the 
smaller region |^i| < A^|^2|- All of these results reduce to the previous case as A — >■ 1. 
In general, these results indicate that the {1}, {2}, and {3} regions become smaller 
as A — > 0, occupying narrower and narrower "pie-slices" in the (^1,^2) landscape 
plane. Specifically, each of these "pie-shces" has total angle 6'{i}^{2},{3} = 2^a, where 

= tan"^ A^ , (3.9) 

and differ only in their orientations in the (^1, ^2) plane, with the Class {1} pie-slice 
centered around the positive ^i-axis, the Class {2} pie-slice centered around the 
^2 = —^1 > diagonal axis, and the Class {3} pic-slice centered around the negative 
(^2-axis. Thus, as A — > 0, the {1}, {2}, and {3} regions disappear entirely. A sketch 
of these results for A = 1/2 is shown in Fig. 4(a). 




Figure 4: (a) Left figure: Same landscape as in Fig. 3, now sketched for A = 1/2. The {1}, 
{2}, and {3} regions have each become smaller, remaining equal in size, while the {12} 
and {23} regions (as well as the new {13} region) have expanded to fill in the gaps. Just 
as for A = 1, the different classes of vacua occupy non-overlapping regions in this landscape, 
and we have indicated the unstable extrema in parentheses, (b) Right figure: As A ^ 0, 
the {1}, {2}, and {3} regions disappear entirely, leaving a landscape entirely dominated 
by non-overlapping {12}, {23}, and {13} regions. These regions ultimately represent flat 
directions in the A = limit. 

Just as in the A = 1 case, extrema in Classes {12} and {23} continue to exist in 
the gaps between the {1} and {2} regions, and {2} and {3} regions, respectively. 
As evident from Fig. 4(a), each of these regions has angle 6'{i2},{23} = 37r/4 — 29 x. 
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Moreover, these extrema are stable everywhere within these regions. Thus, as A — > 0, 
these regions become larger rather than smaller. The new feature for A < 1, however, 
is the emergence of new Class {13} extrema which populate the gap that has opened 
up between the {1} and {3} regions, with angle ^{13} = vr/2 — 26x- Just as with 
the extrema in Classes {12} and {23}, extrema in Class {13} are stable wherever 
they exist. There continue to be no extrema, stable or unstable, in Class {123} for 
A > 0. 

We see, then, that just as for A = 1, the different classes of vacua continue 
to occupy non-overlapping regions for A < 1. Indeed, as A — * 0, the {1}, {2}, 
and {3} regions disappear entirely, leaving a landscape entirely dominated by non- 
overlapping {12}, {13}, and {23} regions, with angles 6'{i2},{23} = 37r/4 and ^{13} = 
7r/2 respectively. This situation is illustrated in Fig. 4(b). 

The A = limit is somewhat special. In this limit, of course, we find that the three 
regions in Fig. 4(b) continue to exist. However, rather than representing stable vacua, 
each of these regions now represents a class of flat directions. This makes sense, since 
SUSY is unbroken for these extrema in the A = limit, even for arbitrary values of 
the FI coefficients ^j. Likewise, it!-symmetry is preserved in these regions. Of course, 
the unstable extrema listed in Fig. 4(b) continue to represent solutions with broken 
SUSY and positive vacuum energies, even when A = 0. 

It is easy to interpret this flat direction in terms of the vacuum expectation values 
Vi. Recall that for any A > 0, there are no extrema in the {123} class. However, for 
A = 0, a new {123} extremum appears, corresponding to the general solution 

l^^ir = ^2!' + ^! , ksl' = \v2?-^2 . (3.10) 

These equations have solutions for (vi, V2, v^) for all (^1, ^2), and thus can be consid- 
ered to exist throughout the entire (^1, C2) landscape. They are also SUSY- preserving, 
maintaining K = for all choices of |f2p in Eq- (3.10), and serve as the algebraic 
representation of the flat direction when A = 0. Indeed, depending on the value of 
(^1)^2), the corresponding {12}, {13}, or {23} solution that appears in Fig. 4(b) is 
only a special case of the more general {123} solution in Eq. (3.10). For example, 
the {12} region in Fig. 4(b), which corresponds to solutions with = .^1 + ^2, 
\v2\^ = and l^ap = 0, can be viewed as a special case of Eq. (3.10) with |f2p = ^2- 
Thus, when A = 0, we should more correctly state that the entire {^1,^2) landscape 
can be described as a single flat direction corresponding to the {123} solution in 
Eq. (3.10). 

3.3 A > 1 

Let us now turn to the situation when A > 1. In such cases, the above algebraic 
conditions for the existence and stability of the {0}, {1}, {2}, and {3} extrema do 
not change. Thus, as A grows beyond 1, we see that the pie-slices occupied by stable 
extrema in these classes actually grow rather than shrink, with angles still given by 
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^{i},{2},{3} = 2^A — 2tan~^ A^. Indeed, as A — > oo, each of these stable subregions 
expands to fill its corresponding existence region, so that all extrema in these classes 
become stable as A — > oo. 

Given these results, we see that Regions {1} and {3} actually overlap for A > 1, 
with overlap angle 6'{i}+{3} = 29x — n/2. Likewise, the {2} region separately overlaps 
with the {1} region and with the {3} region for A > A*, where 

A, = {1 + V2y/^ ^ 1.55 ^ 9, = tan-^A^ = ^ . (3.11) 

o 

In the latter case, these overlap regions have angles ^{i}+{2} = d{2}+{s} = 26'a— 37r/4 = 
2{0x—0^). Thus, we see that we now have cases in which a single point in the landscape 
can give rise to multiple (meta)stable extrema. Of course, only the extremum with 
the lowest energy in each case is truly stable within the specified model, while the 
remaining extremum is only mctastablc. In this example, we find from Eq. (3.7) that 
the {1} vacuum is truly stable for ^2 < Ci; while the {3} vacuum is truly stable for 
< ^2- This conclusion holds for all A > 1. 

We stress that this metastability exists within a single point in the landscape (i.e., 
within a single model characterized by specific values of ^1 and ^2)- This is quite 
different from instabilities across different models at different landscape locations; 
such theory-to-theory transitions are beyond the framework of quantum field theory, 
but have been conjectured to occur in string theory. If such transitions can occur in 
string theory, then we can expect the full dynamics on such a landscape to involve 
both intra-model as well as inter-model metastability. However, there is no evidence 
that such inter-model transitions can occur. 

An important question is to determine the fate of the {12}, {13}, and {23} 
extrema for A > 1. For 1 < A < A*, extrema in the {12} and {23} classes continue to 
exist within their remaining gap regions and are always stable, with angles 6'{i2},{23} = 
37r/4 — 29x. Moreover, for all A > 1, {13} extrema actually re-emerge within the 
Class {l}+{3} overlap region, but are now unstable. As always in this toy model, 
there are no extrema in Class {123} except at A = 0. We thus obtain a landscape 
as sketched in Fig. 5. 

The final remaining case corresponds to A > A*. For A in this range, the {1}, 
{2}, and {3} regions all mutually overlap in a pairwise fashion. (Note that these 
three regions never have a simultaneous common overlap.) In such cases, we then 
find that extrema in Classes {12}, {13}, and {23} all continue to exist in their 
respective overlap regions, but that all of these extrema are now unstable. We thus 
obtain a landscape as sketched in Fig. 6(a). Note that as A — cxd, our landscape 
reduces to that sketched in Fig. 6(b). 



13 



^2 




Figure 5: Same landscape as in previous figures, now sketched for A = 4/3. Regions {1}, 
{2}, and {3} have each become larger, with a new overlap region (indicated with double 
cross-hatching) between the {1} and {3} regions. Class {13} extrema exist but are now 
unstable in this overlap region, while stable Class {12} and {23} vacua continue to popu- 
late the gaps between Regions {1} and {2}, and Regions {2} and {3}, respectively. There 
continue to be no overlaps amongst Regions {l}/{12}/{2} or Regions {2}/{23}/{3}. 
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Figure 6: Same landscape as in previous figures, now sketched (a) for A = 2 (left), and (b) 
as A ^ cx) (right). Since A > A* in each case, the {1}, {2}, and {3} regions all mutually 
overlap; these three overlap regions are indicated with double cross-hatching. Extrema in 
Classes {12}, {13}, and {23} exist in these overlap regions, but are always unstable for 
A > A*. 14 



3.4 Landscape probabilities and integration measures 



Just as in the A = 1 case, we can calculate the fraction of stable extrema relative 
to all extrema as a function of A; we merely evaluate this fraction in each relevant 
pie-slice, and then properly weight these pie-slices according to their opening angles. 
It turns out that as a function of A, this fraction is given by 

' (77r + 10 tan-i A^) /247r for A < 1 , 
Prob (stable; A) = I (IStt + IGtan-^ X'^)/A8n for 1 < A < A* , (3.12) 

, (tt + 2 tan-i A^) /47r for A > A* . 

Note that this probabihty agrees with Eq. (3.6) when A = 1, and generally increases 
monotonically from 7/24 ^ 0.29 at A = to exactly 1/2 at A — oo. For example, 
it is straightforward to verify from Fig. 6(b) that this fraction is indeed 1/2 in the 
A — oo limit. A plot of this probability as a function of A is shown in Fig. 7. Once 
again, we stress that we are averaging over our (^1,^2) landscape with a uniform 
measure when calculating this probability. As already noted below Eq. (3.6), we see 
that this probability is significantly greater for all values of A than we might have 
naively expected. 




Figure 7: The probability that a given extremum in our toy model is stable, plotted as a 
function of the Yukawa coupling A and integrated over all (^1,^2)- This function increases 
monotonically from 7/24 when A = 0, and asymptotes to 1/2 as A ^ 00. Note that in all 
cases, this probability is significantly greater than we would have naively expected based 
on a random assignment of signs for the eigenvalues of the six-dimensional mass matrix. 

Of course, A is actually a landscape coordinate; just as with ^1 and ^2, the value 
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of the parameter A specifies a particular model under study.* In other words, as dis- 
cussed at the beginning of this section, our toy model really has a ^/iree-dimensional 
landscape, with these figures showing only particular fixed-A slices through this land- 
scape. Thus, in calculating the probability of finding a stable extremum from amongst 
all extrema, we should really be averaging over A (our third direction) as well. 

However, this illustrates another critical issue for the landscape. What measure 
should we use when performing our A-integration? If we were to perform our averag- 
ing process using the measure dX, we would find that the resulting average probability 
would be close to 1/2 because of the overwhelming dominance of the contributions 
from large |A|. This suggests that the simple measure dX may not be the most ap- 
propriate measure for the Yukawa coupling A, even with an appropriate perturbative 
cutoff. As briefiy mentioned in Sect. 2, another option prompted by the geometry 
of our landscape might be to consider an alternate measure based on our angular 
Yukawa variable 9\ = tan~^ A^: 

dOx = . (3.13) 

^ 1 + A4 ^ ' 

Note that this measure suppresses the contributions from large |A|, as wc desired, 
yielding a finite result. For example, using this angular measure, we find that the 
average probability for finding stable extrema becomes 

2 /"Vs 151 
Prob (stable) = - / Prob (stable; A) = — 0.393. (3.14) 

71 Jo 384 

One could also implement a perturbative cutoff for 9x, again yielding cutoff-sensitive 
results. However, regardless of the particular measure or cutoffs chosen, this example 
illustrates in a graphic and concrete way the need to have an unambiguous way of 
determining appropriate measures and cutoffs in order to meaningfully calculate any 
statistical quantities on the landscape. 

Another quantity that we can calculate in our toy model is the probability that 
a given (stable) vacuum state preserves it!-symmetry. In our toy model, such R- 
symmetry-preserving vacua are those in Classes {1}, {2}, and {3}, while the vacua 
in Classes {12}, {13}, and {23} have the potential to break ii!-symmetry. Thus, cal- 
culating the probability of randomly obtaining an i?-symmetric vacuum from amongst 
all stable vacua as a function of A, we obtain 

(3/7r) tan-i A^ for A < 1 , 

(1/4) + (2/7r) tan-i A^ for 1 < A < A, , (3.15) 
1 for A > A* . 

This probability is plotted in Fig. 8 function of A. 

■tin a realistic setting, all three of these quantities may change dynamically under renormalization- 
group flow. This issue will be discussed further in Sect. 3.5. However, for the purposes of studying 
the vacuum structure of our toy model, we are considering these quantities as primordial fixed 
parameters defining our specific model. 



Prob (R symm; A) = < 
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Figure 8: The probability that a given vacuum in our toy model preserves i?-symmetry, 
plotted as a function of the Yukawa coupling A and integrated over all (^i, ^2)- This function 
increases monotonically from zero when A = 0, and saturates at 1 for A > A*. 

Once again, just as with the probability of finding stable extrema, we recognize 
that we should really average this result over A (the third landscape coordinate) . If we 
were to use the measure /q°° dX or dX, we would obtain the result that the average 
probability is one, since the contributions from large |A| (or from |A| > A*) dominate 
over the contributions from smaller A. However, using the "angular" measure in 
Eq. (3.13), we obtain the exact analytical result 

2 W2 21 
Prob (R symm) = — / (i6'A Prob (R symm; A) = — ~ 0.656. (3.16) 
TT Jo 32 

Clearly, this result depends strongly on the measure employed, and possible cutoffs 
for A that might be imposed. 

3.5 Renormalization-group flow and boundary crossings 

One of the interesting features of the vacuum structure we are seeing in this toy 
model is the fact that all of the boundaries between different vacuum regions are actu- 
ally energy-dependent (or temperature-dependent in the early universe). Therefore, 
it is possible that a vacuum can cross a boundary between regions as the result of 
renormalization group evolution, either because the landscape location of the vacuum 
changes, because the boundary changes, or both. 

In order to understand this potentially interesting observation, let us now analyze 
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the evolution of our toy model under renormalization-group (RG) flow. As we have 
seen, this model actually contains several quantities which are potentially renormal- 
ized: these include the FI coefficients and the Yukawa coupling A. However, we 
must remember that our toy model also depends on the U{l)i gauge couplings gi that 
were implicitly dropped from Eq. (3.2). In this model, the RGEs for the two gauge 
couplings gi, g2, and the Yukawa coupling A are given by 



where Tr Q"^ is the sum of (squared) charges under U{l)i. 

In order to illustrate the phenomenon of boundary crossing, we will assume for 
simphcity that the entire matter content of our theory consists of the three chiral 
superfields $i, $2, and $3, and that the FI terms (.^1, .^2) are introduced at tree level. 
We therefore find that TiQf = TrQ"^ = 2. Moreover, since TrQi = TrQ2 = in 
this case, the one-loop induced FI coefficients are zero, and therefore the tree-level 
FI coefficients (^1, ^2) are not renormalized. Thus, we see that in this toy model, the 
location of an individual vacuum in the landscape is invariant, and the only changes 
that can occur are those which change the topography of the landscape itself (moving 
boundaries, growing or shrinking regions, etc.) with respect to that fixed location. 

Let us consider initial boundary values 5'^(A) = g^^ and A^(A) = Aq, where A 
is some initial reference ultraviolet (UV) scale. In order to make this calculation 
tractable, let us further assume that our gauge couplings are originally equal at the 
UV scale: gi^o — 92,0 = go- We therefore find from Eq. (3.17) that our gauge couplings 
remain equal at all subsequent energy scales: gi{fi) = g2{.lj) = g{lj)- 

The question that we wish to pursue, then, is that of determining the evolution 
of the entire landscape structure as a function of the energy scale ji. We shall do this 
in several steps. 

First, since g\{[j) — 5'2(a*) = 5'(a*), we see that we can easily restore the gauge 
couplings to all of our previous landscape calculations simply by rescaling A(/i) 

g{li). This holds because the vacuum structure of the potential V in Eq. (3.2) is 
not affected if we balance the rescaling of the D-terms with a corresponding rescaling 
of the F-terms. Of course, such a universal rescahng of A would not be possible if we 
did not assume that gi^ = g2fi- 

Next, let us consider the sequence of landscape "snapshots" that we have already 
generated in Figs. 3, 4(a), 4(b), 5, 6(a), and 6(b). Reordering these figures, we 
see that they describe the landscape as a function of the parameter A(/x), which 
should now be interpreted as the rescaled parameter \{^)/g{^). Therefore, if we 
wish to understand the RG flow in the landscape, we simply need to determine 
the flow of the single quantity \/g which parametrizes which landscape sketch is 
appropriate at which energy scale. Our original location on the landscape doesn't 




(3A2 - Agl - Agl) , 



(3.17) 
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change with energy, but we simply have to look at the correct figure corresponding 
to the appropriate value of \/ g in order to know the physics at this location. This 
observation is important, because it confirms that the RG fiow simply corresponds 
to moving between landscape pictures while holding our original landscape location 
(Ci,o, ^2,0) fixed. There is therefore no need to perform any RG calculation for the 
behavior of each individual boundary region. 

Note that this last observation would not hold if we did not assume that gi^ — g'2,o- 
In the case of unequal gauge couplings, our landscape pictures in previous sections 
might not apply, and new regions of stability/instability might emerge. However, 
because of the universal scaling argument which is possible when (?! = 5'2,0) we see 
that our previous landscape calculations are guaranteed to apply regardless of 

Thus, in order to determine the evolution between different landscape pictures, we 
simply need to determine the RG equation for our rcscalcd Yukawa coupling Y{ii) = 
\{fi)/ g{fi). For example, when gi{n) = g2{fJ') = oifJ'), the boundaries demarcating 
the {1} + {3} overlap region for > 1 can be specified directly in terms of Y{iJ,): 

-Y\^^)i, < 6 < -y^iy^i • (3-18) 

Given Eq. (3.17), it is relatively straightforward to determine the RG equation for 
Y{n). We find 

I^^Y{^,) = -|l^(3y^- 10) /(/.), (3.19) 

ail lOTT^ 

and substituting the explicit solution for 5'^(/x) from Eq. (3.17), we can integrate 
Eq. (3.19) to obtain the solution 
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y2(^) 




9a . 



(3.20) 



Thus, we see that quantity 3 — 10/y^ scales according to the ratio of gauge couphngs 
[9{^A / ■ Since this ratio always decreases as we fiow towards the infrared, we see 
that 3 — 10/y^ always evolves towards zero, either from above or below, as we fiow 
towards the infrared. In other words, regardless of the initial value Yq, our theory 
always fiows towards an infrared fixed point 

— flO 

Y = J— f« 1.826 . (3.21) 

The physical interpretation of this behavior is as follows. We begin our evolution 
with some initial values (which determine our location in the landscape) as well as 
initial values Aq and go (which determine which landscape sketch is appropriate). Our 
landscape location is invariant as we flow towards the infrared, but the underlying 
landscape picture (boundaries, etc.) will evolve. If we find that we are originally in 
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a landscape with Yq ^ \q/ qq < Y , we flow upwards through the landscape diagrams 
towards the landscape with Y — Y . Conversely, if we are originally in a landscape 
with YQ>Y,vfe flow downwards through the landscape diagrams towards the land- 
scape with Y = Y . In all cases, however, the net result of the RG evolution towards 
the infrared is to wind up in the landscape corresponding to Y = Y . This is the 
universal infrared limit. 

Given this, it is straightforward to trace the RG dynamics that a single point 
in the landscape might experience. Note that a generic feature of the boundaries 
separating two different vacuum stability regions is the presence of a masslcss scalar 
in the spectrum. This is intuitively clear, since boundaries demarcate the stability 
regions of various vacua. Thus, crossing a boundary is in some sense equivalent to a 
phase transition, with the appearance of long-range order as the new massless state 
appears in the spectrum. The order of the phase transition can then be determined 
from the behavior of the vacimm energy V as the boundary is crossed. 

Note that in general, there are only two generic classes of boundaries that appear 
in our landscape diagrams for this toy model: 

• Boundaries separating single- vev regions, such as {1}, and two-vev regions, 
such as {12}, in which one of the two vev's is the same as that in the single- 
vev region. These sorts of boundaries appear in Fig. 4(a), for example, and 
tend to dominate at small Y . 

• Boundaries separating single- vev regions, such as {1}, and overlap regions, 
such as {l}-|-{2}, in which one of the overlapping vacua is the same as the 
vacuum in the single- vev region. These sorts of boundaries appear in Fig 6(a), 
for example, and tend to dominate at large Y . Note that near the boundary, 
it is the common vacuum that has the lower vacuum energy V in the overlap 
region, while the other vacuum in the overlap region is only metastable. 

It is easy to verify that the vacuum energy V is continuous across the first class of 
boundaries, while for the second class, the stable and the metastable vacua are not 
degenerate in energy on the boundary. Note that the only other types of boundaries 
are those that separate single-vev regions from each other (such as the boundary 
separating the {1} and {3} regions in Fig. 3), but these are flne-tuned and occur 
only at the critical values Y — 1 and Y — Y^ ^ l.hh [where Y^ corresponds to A* in 
Eq. (3.11)]. 

Given these classes of boundaries, let us therefore consider the kinds of phase 
transitions which can result as a consequence of RG flow in our toy model. 

If Yq > y, then we flnd ourselves in a landscape containing only single-vev and 
overlap regions. Thus, we have only boundaries of the second type. Moreover, since 
RG flow pushes us towards landscapes with smaller values of F, we can only have 
situations in which our overlap regions are getting smaller rather than larger. Thus, 
depending on our original (fixed) landscape location, the only type of boundary 
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crossing that may occur in this case is one in which our location changes from being 
within an overlap region to within a single-vev region. For example, we may be 
originally be located in the {l}+{3} overlap region, closer to the {1} region than 
to the {3} region. In such the boundary between the {l}+{3} region 

and the {1} region moves towards us and passes our location, we ultimately find 
ourselves in the {1} region. 

In this case, the resulting physics depends on our initial vacuum state. Recall 
that in the overlap region, the true stable vacuum is determined according to the 
nearer single-vev region. For example, in {l}+{3} overlap region, the {1} vacimm 
has lower energy than the (metastable) {3} vacuum if we are closer to the {1} 
region than to the {3} region. Thus, if we are originally in the {1} state, then there 
is no effect as the boundary passes our location. However, if we are originally in 
the metastable {3} vacuum, this vacuum becomes unstable and the {1} solution 
becomes the only stable vacuum as the boundary passes our location. Of course, a 
key ingredient in making this determination is the timescale for boundary crossing 
compared with the timescale for decaying from the metastable {3} vacuum to the 
truly stable {1} vacuum even without a boundary crossing. If the latter timescale 
is shorter than the former, the decay from the metastable {3} vacuum to the truly 
stable {1} vacuum can occur even before the boundary reaches us. 

The analysis is slightly more complicated for Fq < ^- If < 1) then our original 
landscape has only one-vev and two-vev regions, with the two-vev regions shrinking 
as a result of the RG flow towards larger F-values. If we are originally located in 
a two-vev region of this landscape, then we will necessarily eventually experience a 
second-order phase transition into a one-vev vacuum as a result of RG flow. However, 
if 1 < Fq < ^* [where corresponds to A* in Eq. (3.11)], then our landscape 
consists of a mixture of one-vev, two-vev, and overlap regions. As a result of RG flow 
towards greater F-values, the two-vev regions are shrinking and the overlap regions 
are growing. Thus, two different types of transitions are possible: either we can be 
located in a two-vev region and experience a second-order phase transition into a 
one-vev region, as described above, or we can be originally located in a single-vev 
region next to a growing overlap region. In the latter case, our vacuum state in 
the single-vev region continues to be the truly stable vacuum state in the overlap 
region, so there is no phase transition. Finally, if <Yq < Y, we flnd ourselves in 
a landscape in which there are only single-vev regions and overlap regions, with the 
overlap regions growing as a result of RG flow. In such a case, as above, no phase 
transitions are possible: if we pass from a single-vev region into an overlap region, 
our original vacuum state continues to be the truly stable vacuum state in the overlap 
region, and no phase transition occurs. 

Clearly, the possibility of such RG-induced phase transitions in the landscape is of 
great interest and relevance for cosmology, and in particular for inflationary models. 

One important consequence of the infrared flxed-point behavior towards Y — 
Y 1.826 follows directly from Eq. (3.15): our theory always flows in the infrared 
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to one in which R-symmetry is preserved. This observation is true in our toy model 
regardless of the original landscape location or Yukawa/gauge couplings. 

Indeed, the emergence of an infrared fixed point has an even more significant 
consequence: in such cases, the low-energy phenomenology becomes insensitive to the 
plethora of (ultimately string-theoretic) variables that define the ultraviolet landscape 
physics. If such infrared fixed points are generic features of the string-theoretic 
landscapes, their existence suggests that it may not be necessary to understand the 
full ultraviolet string theory in order to extract physically testable predictions from 
the landscape. 

Of course, we have made a number of assumptions in performing our analysis in 
this section. In a more realistic situation, or with (71 o 7^ g2,o, the true landscape 
topography is likely to be much more complicated, with an even richer set of possible 
phase transitions. Other important effects that we have ignored include, for example, 
the possibility of kinetic mixing [11]. Similarly, the RG evolution of the gauge cou- 
plings must clearly stop below the scale of gauge symmetry breaking; likewise, the 
beta-function coefficients depend on the scale of supersymmetry breaking in the sense 
that the matter spectrum is supersymmetric above this scale and non-supersymmetric 
below it. Such issues, while very important, are clearly model-dependent and must 
be studied case by case. 

However, even in this simple toy model, we see that RG fiow has great potential to 
lead to a number of novel phase transitions and boundary crossings in the landscape 
— effects which are highly sensitive to our initial landscape location (^1,0,^2,0)- Of 
course, as already emphased, we stress that the possible vacuum-structure phase 
transitions that we have discussed in this section must be understood in terms of 
temperature phase transitions in the early universe. 

3.6 Adding soft mcisses for the chiral fields 

Finally, before concluding this section, let us briefiy comment on the effects which 
are induced in our toy model when non-zero masses are introduced for our chiral scalar 
fields. Towards this end, let us consider adding the mass term 
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to the scalar potential given in Eq. (3.2). 

It is straightforward to repeat the above vacuum-structure analysis for this case 
as well. Let us first focus on the case with A = 0. We then find that Class {0} 
solutions, which were formerly unstable everywhere, now become stable in a small, 
central, triangular region given by 



Likewise, the solutions in Classes {1}, {2}, and {3}, which were also unstable 
everywhere for A = 0, now become stable within rectangular semi-infinite "strips" 




(3.22) 



i=l 



(3.23) 
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emanating outward from this central triangle. The final landscape is shown in Fig. 9. 
Note that as (^i, ^2) ioo {i.e., as we move increasingly far from the origin or either 
axis), this diagram reduces back to that shown in Fig. 4(b). This is in accordance with 
our expectation that the effects of the mass terms should vanish if these masses are 
small compared with the Fayet-Iliopoulos (FI) coefficients. Of course, these masses 
remain important if either FI coefficient is small, or if their sum ^1 + ^2 is small 
with ^1 < 0, 5 > 0. 



^2 




Figure 9: Same A = landscape as shown in Fig. 4(b), but now with non-zero soft masses 
> for the chiral superfields. The non-zero masses open up a central triangular region 
of stable vacua in Class {0}, and likewise introduce three semi-infinite rectangular "strips" 
of stable vacua in Classes {1}, {2}, and {3}. Note that very far from the origin, the 
relative effects of the mass terms generally vanish. In this limit, only those regions with 
non-zero opening angles survive, and the diagram reduces back to that shown in Fig. 4(b). 

For A > 0, the situation generalizes as we might expect. First, we observe that 
stable {0} vacua continue to exist within the triangular region in Eq. (3.23) for all 
A, since this result is A-independent. In general, for any A > 0, the one-vev stability 
regions then emanate outwards from the three sides of this central triangle, with the 
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same opening angles as they have for — 0. Similarly, the stability regions that 
emanate from the vertices of this central triangle are the same as they would have 
been for fh^ = 0, corresponding cither to two-vev vacua or to overlapping one-vev 
vacua depending on the value of A. In all cases, we see that these landscapes with 
> reduce back to the previous landscapes with = when we look far from 
the origin and therefore retain only those stability regions with non-zero opening 
angles. 

Finally, let us briefly consider the case with < 0. In this case, and with 
A = 0, it is easy to verify that no stable vacua exist anywhere in the landscape, i.e., 
for any values of ('^i,^2)- This makes sense, since we already saw that the A = 
landscape comprises a single flat direction when — 0, and this flat direction is 
necessarily destabilized for < 0. However, for non-zero A, this expectation can 
change since we can balance the stabilizing F-terms against the destabilizing soft 
masses. In general, we find that as A increases from zero, stable pie-slice regions of 
stable vacua begin to exist at the outer edges of the landscape. These vacua have 
the same opening angles and properties as functions of A as they have in the = 
case, but their origination points are displaced away from the origin, leaving a hollow, 
unstable center. Thus, this is our flrst example of a hole in the landscape — i.e., a 
bounded finite region of the landscape in which there are no stable vacua of any sort. 

As A increases, each of these pie-slice regions begins to originate from a point 
closer and closer to the origin. Together, they eventually flU the hole completely 
when A reaches a critical value related to |m|. However, as A continues to increase 
beyond this critical value, each of these origination points continues to move past 
the origin, thereby creating a crowded new region of three overlapping single-vev 
vacua in the vicinity of the origin. This is thus our first example of more than two 
simultaneously overlapping vacua in the same landscape region. 

4 Adding more U(l) gauge groups 

In this section we shall consider generalizing the model in Sect. 3 by adding more 
U{1) gauge groups. This will significantly increase the number of vacua and the 
complexity of the corresponding landscape. More importantly, since the Wilson-line 
superpotential can in principle contain more fields, we see from dimensional analysis 
that the F-terms will generically be suppressed. Thus, i?-symmetry will be tend to 
be broken only at very low energies and only for relatively few vacua. 

There are, of course, many different ways in which we can introduce further U{1) 
gauge groups. The simplest example containing a large number of vacua is provided 
by considering n copies of the Fayet-Iliopoulos model analyzed in Sect. 2. Such a 
setup is easily realized in string models containing large numbers of Abelian factors, 
corresponding to branes with significant geometric separations in the six-dimensional 
internal space. Such a model would clearly have 3" classes of vacua, and requires 
the various Abelian gauge group factors to be anomalous in the string-theory sense. 



24 



While this is certainly possible, one must bear in mind that in perturbative string 
constructions, the number of Abelian factors is limited by the total rank of the gauge 
group, as determined by Ramond-Ramond (RR) tadpole cancellation conditions. In 
2" of these classes of vacua, the Abelian gauge group is completely broken.* 
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Table 2: U{1) charge assignments for chiral superfields, inspired by "deconstruction" models 
of extra dimensions. 

Things become less straightforward if the various Abehan factors are not decou- 
pled from one another. For this purpose, let us therefore consider a generalization of 
the two-f/(l) model of Sect. 3 to the case of n different U{1) gauge group factors, with 
n+ 1 chiral superfields. Inspired by deconstruction models of extra dimensions [7], we 
shall take our charge assignments to follow the pattern indicated in Table 2. Thus, 
as evident from Table 2, we are only considering bi-fundamental "nearest- neighbor" 
charges; other configurations will be briefly discussed in Sect. 6. Likewise, we shall as- 
sume for simplicity that only the boundary U{1) gauge-group factors, i.e., U{l)i and 
U{l)n, have Fayet-Iliopoulos coefficients and ^„ respectively. Given these charges, 
we can in general write down a Wilson-line supcrpotential of the form^ 

n+1 

W ^ XH^i. (4.1) 

i=l 

Note that unlike the previous case with n = 2, this supcrpotential is not renormaliz- 
able for n > 2. However, we can continue to study the landscape of this model as a 
function of the Fayet-Iliopoulos coefficients (^i, for arbitrary values of the Yukawa 
parameter A. We can also consider the addition of soft scalar masses by adding a 
term 

Koit = rh'Em' (4.2) 

i=l 

to the scalar potential, as in Sect. 3. 

* Models with several Abelian factors have been considered phenomenologically in various papers; 
see, e.g., Ref. [11] and references therein. 

^This implies that the $i superfields have i?-charge 2/{n + 1), giving the F-terms ii-charge 
—2n/{n + 1). Thus, as claimed earlier, i^-term breaking will correspond to ii-symmetry breaking 
in this model. 
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4.1 Special case n = 3 



Let us first consider the case with n = 3. Since the superpotential W is non- 
renormahzable, its contribution to the overall potential is negligible at energy scales 
below the fundamental physics scale. We shall therefore take A = in our analysis. 
However, we shall continue to leave the common scalar soft mass rh^ arbitrary. 

Since there are four chiral supcrficlds when n — 3, there are now sixteen possible 
classes of vacua which can be characterized by the vacuum expectation values (vcv's) 
Vi of the scalar fields 0j. Wc shall continue to employ a notation in which {ijk...} 
denotes tlic class of extrema in which Vi, vj, Vk, are non-zero. Likewise, the 
notation {0} will denote the class of extrema in which all vev's vanish. For example, 
{14} denotes extrema in which V2 = = while Vi 7^ 0, V4 7^ 0. In general, just as 
in the n = 2 case, we can calculate a (8 x 8)-dimensional mass matrix as in Eq. (3.5). 
There will be a zero eigenvalue for each spontaneously broken these eigenvalues 

correspond to the resulting Nambu-Goldstone bosons. Our extrema then represent 
stable (or metastable) vacua only if each of the remaining eigenvalues is positive. 

We find the following results. In general, of the sixteen potential classes of ex- 
trema, we find that only eleven classes of extrema actually exist; there are no solutions 
for extrema in Classes {12}, {34}, {124}, {134}, or {1234}. 

Of the remaining eleven classes of extrema, we find that only three of these give 
rise to stable (non-overlapping) regions when m — 0. This case is sketched in Fig. 10, 
which should be compared with the analogous n — 2 case Fig. 4(b). Compared with 
the n = 2 case, we see that the stable regions have a virtually identical shape, with 
the {12} vacua for n = 2 corresponding to the {123} vacua for n = 3; the {23} 
vacua for n = 2 corresponding to the {234} vacua for n = 3; and the {13} vacua 
for n = 2 corresponding to the {14} vacua for n — 3. There are some differences, 
however. For example, in the n — 2 case, the entire landscape was dominated by 
vacua that break i?-symmetry. By contrast, if we classify our n = 3 vacua according 
to their general properties for arbitrary A, we see that the {14} vacua are actually 
supersymmetric and consequently preserve i?-symmetry. This is therefore our first 
example of a landscape containing both supersymmetric and non-supersymmetric 
regions of the same dimensionality [as opposed to lines of supersymmetric vacua 
within two- (or higher-) dimensional non-supersymmetric regions]. We also notice 
that the n = 3 case gives rise to slightly more unstable extrema for each stable 
vacuum than we found in the n = 2 case. Finally, we note that the dashed lines in 
Fig. 10 have slopes 2 and 1/2, as opposed to the single slope 1 in Fig. 4(b). 

We stress that Fig. 10 must be interpreted in the same way as Fig. 4(b), namely 
as the landscape that emerges for infinitesimal A approaching zero. When A = 0, 
the entire landscape in Fig. 10 actually becomes supersymmetric and i?-symmetry 
preserving, completely dominated by the flat directions implicit in the A = solution 
{1234}. 

We can also consider the effects of turning on the soft masses rh. As expected 
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Figure 10: The landscape of our n = 3 model for m = 0. Stable vacua are indicated in 
large bold type, whereas unstable extrema are indicated in smaller non-bold type. Vacua 
breaking i?-symmetry are shaded; note that the {14} region is supersymmetric and hence 
i?-symmetry preserving. Compared with the analogous n = 2 case in Fig. 4(b), we see that 
not all vacua break i2-symmetry for n = 3. There are also slightly more unstable extrema 
for each stable vacuum. 

from the n — 2 case, we find that the {0} solutions, which were formerly unstable 
everywhere, now become stable in a small, central, square region given by 

161 161 <m'. (4.3) 

Likewise, of the seven remaining unstable classes of extrema, two of these (Classes 
{2} and {3}) become stable within small triangular regions attached to this central 
square, while the remaining five classes of vacua now become stable within rectangular 
semi-infinite strips emanating outward from the central square and triangles. The 
resulting landscape is shown in Fig. 11. Detailed solutions for each class of extrema 
in this model are given in the Appendix. 



27 



^3 




Figure 11: The landscape of our n = 3 model, now drawn for > 0. Stable vacua are 
indicated in large bold type, whereas unstable extrema are indicated in smaller non-bold 
type. Vacua breaking i?-symmetry are shaded. This landscape is to be compared with its 
n = 2 equivalent in Fig. 9. Once again, we observe that in the m — limit {i.e., as we 
move radially outwards from the origin), only those regions with non-zero opening angles 
survive; this diagram then reduces back to that in Fig. 10. 



4.2 Arbitrary n and large-n limit 

We now turn our attention to the general-n case, with particular interest in the 
large-n limit. In the general case, the D- and F-terms are now given by 

Da ^ + il^al + iJan , ^ —- , (4.4) 

where a — 1, . . . ,n. This leads to the scalar potential 

n n+1 

V = lEDl + ■ (4-5) 

a=l i=l 
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As we already remarked for the n — 3 case, the F-term contributions coming from 
the Wilson line superpotential (4.1) in the large-n case have a neghgible effect on 
the vacuum structure for FI terms smaller than the Planck (or string) scale. The 
extrema of the scalar potential are then obtained as solutions to the field equations 



0j 



/ ^1 -.1 \ 

n n+1 n+1 

Oi'^^l + qt^n + E E + I E I<^1 ■ ■ ■ ■ ■ ■ 4- ■ ■ ■ <t>n+l? 



= 



(4.6) 

where the hatted fields refer to fields that are omitted from the product. These 
equations each have two solutions: = {4>k) = 0, and a solution obtained by setting 
the term in parentheses to zero. This gives a large number ~ 0(2"^+^) of different 
extrema which were explicitly analyzed in Sect. 3 for n = 2 and Sect. 4.1 for n = 3. 

Certain features of the general-n solution can easily be perceived directly by gener- 
alizing from the explicit n = 2 and n = 3 solutions. For example, it is straightforward 
to show that the m = 0, A = landscape always has three non-overlapping regions 
with the same shapes as in Fig. 6(b) and Fig. 10, and with the dashed lines (sepa- 
rating regions with different unstable extrema) having slopes n — 1 and l/{n — 1). 
Moreover, it is easy to sec that the three regions of stability always correspond to 
vacua in the classes {l,n-|-l}, {l,2,3,..,n}, and {2,3,...,n-|-l}. 

Similar general statements can also be made for the m 7^ case. When > 0, 
the presence of the non-zero soft mass always opens up a central region of stability 
for the {0} extrema. In general, this region is an (n-l-l)-gon, with regions of stability 
for otherwise unstable extrema either emanating outwards in semi-infinite strips, or 
attached to the {n + l)-gon with bounded regions shaped like n-gons. Note that the 
appearance of a central (n + l)-gon can be understood as the result of having n + 1 
independent eigenvalue constraints defining the borders of the stability region. 

While all of these properties hold true for general n, it is also instructive to exam- 
ine certain specific vacuum solutions. Instead of attempting a general classification, 
we will concentrate on presenting some general features such as the absolute minima 
of the theory. 

In the following, we shall focus on the case in which \/M'^-^ ^ and m = 0. 
Since X/Mp""^ — > 0, we expect that the stable vacua are dominated by the three 
classes {l,n-|-l}, {l,...,n}, and {2,...,n-|-l}. Specifically, we find that the {l,n-|-l} 
solution exists for ^1 > 0, ^„ < 0, and is given by 

{l,n + l}: |-f;i|2 = ^^, |^^^^|2^_^^^ ^2^___^^^^0. (4.7) 
Similarly, the {l,...,n} solution exists for ^„ > 0, ^1 -|- ^„ > 0, and is given by 
{l,...,n}: KP = (6 + en)(l + 0(6^)) 

\V2\\ M = en(l + 0{el)) , Vn+1 = (4.8) 
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where ~ (A/M^"^)2(^„/M|,)"-2, while the {2,...,n+l} solution exists for < 0, 
+ < 0) is given by 

{2, n + 1} : v, = 0, \vn+i\' = -(6 + Ui^ + 0{ei)) 

\v2\',...,\vn\' = -U^ + 0{ei)) (4.9) 

where ej ~ {\/M'^-y{Ci/Mj,)''-'^. Note that ei and e„ are very small numbers 
for n ^ 1 and C,i ^ Mp. The vacuum in Eq. (4.8) was discussed in detail in 
Ref. [12]. Each of these classes of vacua occupy non-overlapping regions in the two- 
dimensional (C,Cn) parameter space. Note that the {l,2,...,n} and {2,3,...,n-|-l} 
vacua completely break all of the U{1) gauge factors, while the {l,n-|-l} vacuum is 
supcrsymmetric and A-independent for all n > 3. 

We shall also be interested in several other explicit solutions for general n. All of 
the following solutions are A-independent. For example, the {2,3,...,n} solution is 
given by 

{2,3, ...,n}: vi ^ Vn+i ^ , 

\vk\''^-[{^i + Q{k-l)-n^i] for A; = 2,..., n. (4.10) 

This solution has an unbroken U{1) generator Qi + ... + Q^, where Qi are the gen- 
erators of U{l)i, and gives rise to the D-terms 

{D,) = ... = {D^) = ^^±^ . (4.11) 

Note that this solution has a linear "profile" , in the sense that the sequence of non- 
zero vacuum expectation values in Eq. (4.10) grows linearly with k. 
Another solution of interest is the {2,3,...,n-l,n-|-l} solution, given by 

{2,3, ...,n- l,n + l} : vi ^ ^ , {vn+il"^ ^ -^n 

\vkf = ^^^i ior k = 2,...,n-l . (4.12) 
n — 1 

This solution has the unbroken generator Qi + Q2 H 1- Qn-i, and gives rise to the 

D-terms 

(D,) = ... = (L>„_i) = -A_ , ^ . (4.13) 

This solution also exhibits a hnear profile. 

Another solution of interest is the {1,3,4, ...,n} solution, given by 

{1,3,4, ...,n} : V2 = v^+i = , \vi\'^ = 

k — 2 

\vk?^ Tin for A; = 3,4,...,n . (4.14) 

n — 1 
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This solution has an unbroken U{1) generator Q2 + ■■■ + Qn, and gives rise to the 
D-terms 

(A) = 0, (D,) = ... = (L>„) = (4.15) 

Clearly, this solution also exhibits a linear profile. 

Finally, the last solutions we shall consider are the {2,3,...,j-l,j+l,...,n} solu- 
tions, given by 

{2,3, - l,j + 1, ...,n} : vi = Vj = Vn+i = 
I^Ci for^ = 2,...,i-l 



n-j + l ^^ for A;=j + l,...,n . 



(4.16) 



These solutions leave unbroken the two generators Qi + Q2 + ■■■ + Qj-i and Qj + 
Qj+i + ... + Qn, and gives rise to the D-terms 

(Di) = ... = = Ci/O'-l) = d, 

{Dj) = ... = (L>„) = ^„/(n-j + l) = (4.17) 

Unhke the previous solutions, each of these solutions has two independent hnear 
profiles. 

Needless to say, there are numerous other solutions which can be generated for 
general n. However, the above solutions will be sufficient for our purposes. 

Note that these vacua, as well as all vacua with smaller numbers of non-zero vev's, 
are unstable at the level of our discussion (in which we are taking A/M^-^ and 
rh — 0). Moreover, one of the obvious properties of such solutions is the fact that, 
at first sight, they unfortunately appear to give rise to a large number of unbroken 
Abelian gauge symmetries which would survive down to low energies. However, as 
we shall see in Sect. 4.4, this is not the case due to the presence of mixed gauge 
anomalies. 

4.3 Higher-dimensional flux interpretations 

We shall now demonstrate that many of the above general-n solutions have natural 
interpretations in terms of higher-dimensional fiux compactifications. 

First, let us consider the {l,2,...,n} and {2,3,...,n-|-l} vacua. These clearly 
can be interpreted as emerging from a five-dimensional supersymmetric U{1) gauge 
theory compactified on the orbifold S^/'^2, with compactification radii 



R ~ n/^^n and ~ n/Vl6| , (4.18) 

respectively. In each case, the four-dimensional zero-mode gauge field receives a 
mass from the four- dimensional FI term ^1 + ^n- The supersymmetry-breaking scale 
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in these two cases is controlled by the Wilson-hne supcrpotential, and all soft masses 
are ~ 0{en) in the first case and ~ 0{€i) in the second case. By contrast, the 
third vacuum {l,n + 1} is supersymmetric and probably does not have an extra- 
dimensional interpretation. 

We now turn to the {2,3,.. .,n} solution in Eq. (4.10). We shall now argue that 
this vacuum can be given the higher-dimensional interpretation of having magnetic 
flux on a torus in a six-dimensional Abelian gauge theory. Indeed, as we will see, the 
smoking gun for such a magnetic Atix interpretation is the presence of a linear proflle 
in the associated vacuum expectation values. 

From a five- dimensional IR^ x S^/'2.2 perspective, where the fifth dimension is the 
interval < y < ttR, a supersymmetric Abelian gauge theory contains a gauge field 
and a Z2-odd real scalar E. The D-term from the four-dimensional point of view in 
the continuous limit of the deconstruction setup discussed above is given by 

L> = SsS + {-\(p,\^ + ^,)S{y) + {\(p^^,\^ + ^^)S{y-nR) . (4.19) 

The standard profile for the scalar E, largely discussed in the literature, is of the 
form (S) = e(i/)^„/2, which in the case + = is the needed profile for pre- 
serving supersymmetry and the gauge symmetry [13]. Notice, however, that the field 
equations 

5E |a|E + ^5 [asE + (Ci - \v,\'')5iy) + (e„ + K+il^iy - ttR)] | = (4.20) 
have another solution on the orbifold 5*1 /Z2, namely 

This solution does not describe the absolute minimum of the theory, since it has a 

large positive vacuum energy, but it is an extremum of the theory. By using ~ n, 
it is clear that Eq. (4.21) matches the deconstructed result Eq. (4.10). On the other 
hand, from a six-dimensional perspective, E corresponds to the sixth-component Aq 
of the gauge field. Then the flux in the two-dimensional compact space is given by 

F56 = (95E) = - - CnSiy-nR) . (4.22) 

The first term is the magnetic flux we were searching for, whereas the localized terms, 
already discussed in the literature, have the interpretation of fluxes localized at the 
orbifold flxed points. Note that the integrated flux in the compact space is actually 
zero, 

/ttR 
dy (^sE) = , (4.23) 
-ttR 

the magnetic flux cancelling the localized contributions at the flxed points. 
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The other extrema hsted in Sect. 4.2 also have a flux-compactification interpreta- 
tion in extra dimensions. For example, let us consider the {2,3,...,n-l,n+l} solution 
given in Eq. (4.12). In the extra-dimensional interpretation, the bulk field E has a 
difi^erent profile: 

^ ' 27rR 2 ^ ' ^ ' 1 for y = 7ri? . ^ ' 

This gives rise to a fiux which again integrates to zero: 

= {d,T.) = TTB- my) , r (^sS) = . (4.25) 

It may seem strange that the discontinuity aX y = 7iR in the auxiliary field D 
can be consistent with the field equations (4.20). However, in Eq. (4.24) we have 
(S(7ri?)) = 0, and by choosing the consistent boundary condition 6T,{7iR) = 0, the 
field equations (4.20) are indeed satisfied. 

This will be the general rule for the fiux interpretation of the solutions with lower 
numbers of non-zero vev's: discontinuities in the auxiliary Z)-terms (and therefore 
the flux) correspond to points where the wave function of S has Dirichlet boundary 
conditions. A solution with i zero vev's for the link fields therefore corresponds to a 
wave function of E having i nodes (and corresponding Dirichlet boundary conditions) 
on the corresponding points. The solution in Eq. (4.12) is also unstable, having a 
tachyon at the position (0„) of the link field having a zero vev, but by adding the 
Wilson-line superpotential contributions, this solution becomes a stable vacuum close 
to the vertical axis = 0. 

The final example with n — 1 non-zero vev's is the {1,3,4, ...,n} solution given in 
Eq. (4.14). In the extra-dimensional interpretation, the bulk field E has the profile 

/y\ _ /r,\ _ /W(27ri?) for < y < vri? , . 

- M^' " lO fory = 0, ^^"^^^ 

leading to a fiux 

^56 = ^ - CnS{y-TTR) , dy (^sE) = . (4.27) 

znR J-ttR 

Moreover, consistent with our previous explanations, (E) vanishes at the position 
y — where the auxiliary field is discontinuous. Just as with the two previous fiux 
solutions, this solution is unstable. However, by adding the Wilson-line superpoten- 
tial, this solution becomes stable close to the horizontal axis = 0. 

The last examples we shall discuss are the {2,3,...,j-l,j-|-l,...,n} solutions. In 
these cases, the profile for the bulk scalar E satisfies 

SJ:!^{d5^)+CiS{y) + CnS{y-7rR)^ 

= 5E + 1(^2 -c^i)[e(l/-%)-e(l/ + %•)]}, (4.28) 
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where yj is the location at which the auxihary field {D) jumps from {D) = di for 
< \y\ < yj to (D) — ^2 for yj < \y\ < nR. The profile has the form 



serves as the continuity condition for the profile at y = yj. By matching the decon- 
structed version to the continuous version via yj ~ j — 1, we find the expected result 
that {E{yj)) — 0. The flux in this case is 



^56 = d2 + lid2- di) [e{y - yj) - e{y + %•)] - 6 5{y) - Cn S{y - ttR) , (4.31) 



which again integrates to zero. 

4.4 Anomalies in the deconstructed C/(l)" theory 

As we remarked earlier, the solutions with n — 1 non-zero vcv's — and indeed all 
of the solutions with fewer non-zero vev's — are unstable at the level of our discussion 
(in which we are taking A/Mp"^ — > and rh — 0). However, at first sight it may also 
seem that these vacua have a further unpleasant feature, namely a large number of 
unbroken Abelian gauge symmetries which survive down to low energies. However, 
as we shall now demonstrate, this is not the case thanks to the presence of mixed 
gauge anomalies in such theories. 

It is easy to see why subtleties related to mixed anomalies will be present in 
any model containing multiple Abelian gauge fields related to higher-dimensional 
compactifications and/or string theory. The basic reason is that bi-fundamental 
fields, which have been the building blocks of our toy landscape models, create mixed 
anomalies which typically require axionic fields and four-dimensional Green-Schwarz 
mechanisms for their cancellation. The axionic fields will come in complex (moduli) 
superfields due to supersymmetry, which will also fix the values of the FI terms of 
the model. The dynamics and the stabilization of these moduli is, of course, very 
important in the landscape picture. Indeed, we have thus far been assuming that 
these moduli are already stabilized by the dynamics (through, e.g., NS-NS and RR 
fiuxes, gaugino condensation, etc.). 

In the C/(l)" example, the mixed anomahes are described by the {n x n)- 
dimensional anomaly matrix 




for < \y\ < yj 
for yj < \y\ < ttR , 



(4.29) 



where the constraint 




(4.30) 





(4.32) 
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where Qa is the generator of the U{l)a gauge group. The exphcit matrix elements of 
Cab are given by [12] 



Cab — 



47r2 



Sa,b-1 - Sa,b+1 + T:tQlSa,lSb,l + Ti:QlSa,nSb,n , (4-33) 



where for generahty we have aUowed for additional particles {e.g., MSSM particles) at 
the fixed points. Otherwise, in our minimal example discussed earlier, we would have 
TrQl = TrQ^ = 0. Whereas at the four-dimensional level the anomaly-cancelling 
mechanism is not uniquely fixed (and we shall shortly provide a trivial illustrative 
example of this), the requirement of realizing a five-dimcnsional Lagrangian in the 
continuous limit is much more constraining. This was consistently done in Rcf. [12], 
where it was shown that the axionic-type fields taking care of the mixed anomalies in 
Eq. (4.33) are themselves bi-fundamental. More precisely, there arc axionic couplings 
of the type — ln{(f)i(f)j)SiaSj^a+iW^Wa,a, in superfield notation. There are also two 
boundary axionic fields Sl,r which transform only under the first and the last gauge 
group. * 

The presence of the non-minimal kinetic (and Kahler) functions clearly modifies 
the minimization procedure. It is therefore important to verify that our previous 
results remain intact. We have explicitly verified that the absolute minima of the 
model, as well as the first (n — l)-vcv example with a clean magnetic flux interpre- 
tation, are affected in the full theory only at subleading order in {j/Mp. However, 
the vacua with fewer than {n — 1) non-zero vev's do not remain extrema in the full 
theory, since a vanishing vev for a bi-fundamental represents a singular point of the 
theory. The natural solution in the case where a zero vev occurs in a link flcld (pk is, 
in the extra-dimensional interpretation, to add a brane at the corresponding point in 
the compact space. We should also introduce a superfield Sk with the gauge trans- 
formation 5Sk ~ Afe_i — Afe, and replace In^^ — > Sk in the gauge kinetic functions. 
This new field would then have the interpretation of couphng to the U{1) gauge field 
of the brane. This automatically cancels the anomalies, explains why a new massless 
gauge field appears in the spectrum, and gives a rationale for the Dirichlet boundary 
condition which we found in Sect. 4.3 at the position of the vanishing vev. 

This new axionic coupling can, however, render the new gauge field massive. 
Indeed, more and more vanishing vev's require the addition of more and more U{1) 
branes and corresponding new axionic couplings, which eventually render most of the 
new brane [/(l)'s massive. 

If we completely abandon the five-dimensional interpretation, there is certainly 
more freedom in anomaly cancellation. For example, we can introduce n moduli 
fields Sa which have U{l)a gauge transformations SSi ~ — A2, SSa ~ Aa_i — A^+i 
for a = 2, ...,n — 1, and SSn ~ A„_i where A^ is the superfield describing U{l)a 
gauge transformations. Coupling these moduli fields to the gauge fields through a 
term of the form J2a SaW^Wa.a will indeed cancel the mixed anomalies. In this case, 

•l-Further details and a discussion of the five-dimensional limit are given in Ref. [12]. 



35 



assuming an a priori stabilization of moduli, all the results of our toy model follow 
accordingly. Generically however, this generates FI terms at all sites {i.e., for all 
C/(l)'s), unless the Lagrangian of the moduli fields and their vev's are such that FI 
terms vanish at all sites but the first and last. Indeed, such examples are relatively 
easy to find. However, the price to pay in such a construction is the loss of a five- 
dimensional interpretation in the limit of large n, and the unsolved remaining issue 
of moduli stabilization fixing the FI terms. 

The upshot of the present discussion, then, is that in order to cancel gauge anoma- 
lies, we must introduce new axionic (supcr)ficlds with appropriate couplings imder 
the various gauge groups. In this case, most of the U{1) gauge-group factors become 
massive, a welcome feature for the phenomenological viabihty of these models. 

5 Extension to supergravity 

Among the main motivations for seriously considering the landscape picture of 
string theory arc its ramifications for vacuum stability and the cosmological constant 
problem. From the perspective of both of these problems, the inclusion of supergrav- 
ity interactions is of crucial importance. We shall now briefly discuss each of these 
issues in turn. 

5.1 Vacuum stability, soft masses, and split supersymmetry 

If we consider only the Abelian gauge flelds which are relevant for our purposes, 
the scalar supergravity potential is given in terms of the Kahler potential K and the 
superpotential W by [14] 

= (K'~WiWDjW - 3\W\'^) 

+ IT. fab {dtU'Ki + O (gf Vi^i + a) (5.1) 

where we have temporarily set Mp — 1 and where we have employed the usual 
notation Ki = diK, Kfj = didjK, and Fi = DiW — diW + KiW. Here fab arc the 
gauge kinetic functions of the Abelian gauge groups, which depend on the moduli 
fields which we are implicitly assuming to be stabilized. Of course, after moduli 
stabilization, we should be considering gauged i?-symmetries at the supergravity 
level [15]. However, since our main focus is on the scalar potential, this subtlety does 
not affect the present discussion and it is enough to consider the expectation value 
ifab) = (^ab/Oa- Notc that the gravitino mass is = |M^pexp(i^). 

As before, the stability of a given extremum in the supergravity framework de- 
pends on the scalar mass matrix: an extremum represents a stable vacuum or flat 
direction only if the mass matrix has no negative eigenvalues. Of a particular im- 
portance for our purposes are the D-term contributions to the scalar masses, which 
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after moduli stabilization are given by 



did^Vo = glq^im + ^KIj)Da + Qik(p'^KiK^ 



(5.2) 



where the charge matrix Qij is defined as 



(5.3) 



a=l 



The supergravity effects on the C/(l)" deconstruction model described in the pre- 
vious section can be easily described in the limit in which the supergravity moduli 
fields are heavy and can be integrated out. This leads to a low-energy theory with 
soft breaking terms. It is clear from gauge invariance that the resulting soft breaking 
terms are of the form 



The 74-term induced by the Wilson-line superpotential has a negligible effect on the 
dynamics for large n, and can safely be neglected in the following. However, the 
effects of the soft terms on the vacua described in previous sections are quite 
clear. Specifically, just as we saw in previous sections, the large number of fiux 
vacua which were previously unstable become now stable in a certain region of the 
parameter space (^i,^n)- 

Clearly, this method of increasing the number of stable vacua only works for 
large soft masses, of the order of the Fl terms. From this point of view, a large 
supersymmetry breaking scale stabilizing most of the landscape vacua fits well with 
the "split supersymmetry" proposal in Ref. [16], in which a large supersymmetry- 
breaking scale was assumed to go together with the presence of a large number of 
vacua. 

5.2 A supergravity toy model and the cosmological constant problem 

In order to quantify the issue of the cancellation of the cosmological constant in 
a simpler setting, we shall now propose a simple supergravity model, inspired by the 
toy models discussed in the previous sections and by the fiux compactifications of 
string theory [6]. Since the deconstruction model of Sect. 3 is difficult to analyze in 
this context, we will instead consider the supergravity extension of n copies of the 
original Fl model. Thus, our model will be based on the gauge group U (1)" and chiral 
superfields with charges ±1 under the gauge group factor f/(l)j. These 

fields are therefore coupled to each other only through supergravity interactions. 



n+1 n+1 



(5.4) 
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The model has a single Kahler modulus T describing the volume of a Calabi-Yau or 
orbifold compactification, and has supergravity potentials described by 




W ^ Y^^i + W^o + ^ 6"''°^ , (5.5) 



i=l 



where Wq is determined by the fluxes according to Ref. [6] and the nonperturba- 
tive (exponential) term in the superpotential, responsible for the stabilization of the 
Kahler modulus T, could come from D7-branes, gaugino condensation, or Euclidean 
D3-branes [17]. All other moduli are assumed to be stabihzed by other dynamical 
effects, e.g., by adding appropriate fluxes. 

If we flrst set the Fl terms to zero, i.e., = 0, there is a supersymmetric solution 
giving an AdS space, obtained by solving = 0. However, in the presence of the 
FI terms, the vev's of the moduli flelds will be changed. If the Fl terms and the 
vev's of the flelds (f)j^^ are sufficiently small compared to the Planck scale, we can 
effectively integrate out the T-fleld. This will give rise to explicit soft terms in the 
Lagrangian [18]. The part of the scalar potential which is relevant for discussing the 
cosmological constant in the various vacua of the theory is then given by 



71 2 

+ E l4"¥ + e.)^ - ^ml/M , (5.6) 

i=i ^ 



where and Bi describe the soft terms and where we have assumed equal soft 
masses for the flelds ^j"*"^ and (^^^ for simphcity. 

A few words of caution are in order here. Note that the potential (5.6) has been 
written to the lowest order in the matter flelds 0^^^. Thus, for example, the gravitino 
mass, which is a function of all the flelds, becomes a constant. As a result of these 
approximations, the vacua derived from Eq. (5.6) are certainly only approximations 
to the true vacua. Note, in addition, that we have also neglected perturbative loop 
effects as well as non-perturbative effects. Moreover, integrating out the held T in 
Eq. (5.5) is a good approximation only if the flcld T is heavier than the charged flelds 
0j. This is not always a good approximation to the dynamics of the model. Despite 
these issues, we believe that our discussion concerning the possible cancellation of 
the cosmological constant will still be of qualitative value, and that the qualitative 
features of such a cancellation — which relies on balancing the positive contributions 
to the vacuum energy coming from the super symmetry-breaking against the negative 
contributions coming from the supergravity interactions — capture the qualitative 
features of the full quantum problem. 
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Minimizing the scalar potential in Eq. (5.6), we find 3" possible vacua, determined 
as simultaneous solutions of the coupled equations 

{m^ + m^^+g^Di)vi^^ + BiV^'^ = 

{mj + m^-g^Di)vl-^ + b4+^ = (5.7) 

where ff^^ = {(pi^^)- Indeed, the field equations (5.7) have three solutions for each i. 
By performing U{1) rotations, the solutions can always be made real. 
First, it is useful to define 



m% = m^ + 7h^, A^ = ^m%,-Bf (5.8) 
for notational simplicity. We then find that the first class of solutions are given by 
_ (Af + my(Af + fffeO I (-). _ {Aj - m%){Al + gU,) 

which leads to D-terms and vacuum energy Aq of the form 

gUD.) = - A^ Ao = ^ ^ (^-i - ^) - 3 m^/,M| . (5.10) 
The second class of solutions, by contrast, are given by 



_ i^t - m%m-gUi) I (-)p ^ _(A[+m|J(Af - g^i,) 



I - r,„2 A2 ' \ - o„2 A2 ' y^-^^) 

which lead to 

gm) = A? , Ao = E # i-\ + ^ml/,Ml . (5.12) 

Clearly these solutions exist only when Ifi^"*!^ > 0, which means that g}\ii\ > A^ and 
m| j > A^ where < for Eq. (5.10) and > for Eq. (5.12). One can also verify 
that these solutions are stable. Finally, the third class of solutions is trivially 

vt^ = vt^ - , (5.13) 
which is easily verified to be stable if 



rr4^,±^Bf + gta > 0. (5.14) 

Of course, in the limit of vanishing soft terms, i.e., mf, Bi 0, all three classes of so- 
lutions reduce to the standard global minima of the Fayet-Iliopoulos model discussed 
in Sect. 2. 
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Prom a string-theory perspective, the hmit of interest corresponds to taking ^ 
ml^mf^Bi. In this limit, the cancellation of the cosmological constant then takes the 
approximate form 

n 

^|e.|VK + mf)2-i?2 ^ ^ml/^Ml , (5.15) 

i=l 

for the 2" vacua in the classes in Eqs. (5.10) and (5.12). In gravity-mediated super- 
symmetry breaking scenarios such as we are considering, we expect rnf,Bi ~ '>ti1^2- 
Since typical values for the Fayet-Iliopoulos coefficients are smaller by one or two 
orders of magnitude compared to the Planck scale, we see that there are therefore ba- 
sically only two possible ways in which these 2" vacua can yield a small cosmological 
constant: 

• We can have 3> my 2- However, in this case, partial cancellation can oc- 
cur only for a small number of gauge factors, which means that a reasonable 
accommodation for a small cosmological constant is difficult to achieve. 

• Alternatively, we can have < ml^2- this case, the cosmological constant 
can be cancelled by combining the two different types of vacua in Eqs. (5.10) 
and (5.12). For example, we would require n ~ 10 — 100 for generic values 

~ (10~^ — 10~^)Mp . Indeed, in such cases, there will be a large number of 
landscape vacua with small vacuum energies. Of course, any vacuum solutions 
which also contain the third class of solutions in Eq. (5.13) will add contri- 
butions ^gj^j to the vacuum energy Aq. These contributions are very large, 
and will upset the cancellation of the predicted cosmological constant unless 

161 <^3/2Mp. 

By contrast, in the opposite limit of small FI terms, cancellation occurs for typical 
values of soft terms m^ ~ \/3m3/2Mp/n. 

The outcome of this discussion, then, is that it is the second possibility itemized 
above which can best accommodate a small cosmological constant in our toy model. 
Barring additional fine-tunings, this clearly points towards requiring <^ m^, from 
which we conclude that even a qualitative supergravity analysis of the vacuum struc- 
ture of any landscape model is crucial for addressing the smallness of the cosmological 
constant.* 

High-scale supersymmetry breaking therefore seems to have two positive features 
in our D-term landscape models. First, they render most of the vacua stable by 
adding large diagonal soft masses. Second, they allow a possible cancellation of the 
cosmological constant. 

*0f course, in addition to the stability of the vacua and the fine-tuning of the cosmological 
constant (both of which seem to favor large soft masses), a second crucial fine-tuning in this case 
is the light Higgs mass. This can only be addressed within the context of a specific model in which 
Standard Model or MSSM fields are coupled to the U{1) fields we are investigating, and is thus 
beyond the scope of the current paper. 
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6 Discussion 



In this paper, we proposed a field-theoretic framework giving rise to models con- 
taining large numbers of vacuum solutions. The field-theoretic nature of these models 
therefore allowed us to explicitly calculate such quantities of interest as the ratio of 
stable versus total numbers of vacua, the number of i?-symmetry preserving vacua, 
and the supersymmetry-breaking scale. While at the field-theory level many models 
with large numbers of vacua can certainly be given, our examples have the advantage 
of describing large classes of string compactifications with Abclian gauge groups and 
FI terms. Moreover, within this large class, we presented specific examples involving 
discretized versions of magnetic fluxes in the internal space, as obtained by decon- 
structing (supersymmetric) models with U{1) gauge fields on the orbifold 3^/^.2. By 
examining the extrema involving vanishing vev's for bi-fundamental fields, we found 
that these solutions correspond to profiles which are linear in for the odd-scalar 
E = ylg in the five-dimensional vector multiplet. These solutions therefore corre- 
spond to constant magnetic fluxes d^T^ — F^q from a six-dimensional perspective. 
The number of deconstructed magnetic fluxes is large, and although tlicy all have 
Nielsen-Olesen instabilities, they can be stabilized by F-term contributions coming 
from supcrgravity interactions. This therefore generates a field-theory landscape 
which is similar in spirit and closely related to the landscapes currently under dis- 
cussion in string-theory contexts [2]. We also showed that (super) gravity effects turn 
out to be consistent with (and actually needed for) a qualitative accommodation of 
the smallness of the cosmo logical constant. 

One of the interesting results of the landscape picture emerging in the class of 
models we considered in this paper is the possibility of passing from one vacuum to 
another by renormalization group flow. This possibility arises because our funda- 
mental deflning parameters, such as and A, can change with the energy scale. We 
showed in an explicit toy model that this renormalization group flow can indeed in- 
duce boundary crossings. A general feature of the boundary separating two different 
vacua is the presence of a massless scalar in the spectrum, which suggests an interpre- 
tation in terms of phase transitions. This type of boundary crossing could therefore be 
of potential interest for inflationary models. Moreover, the gauge symmetry-breaking 
pattern can also change as a result of boundary crossings. The corresponding phase 
transitions could therefore also be relevant for electroweak symmetry-breaking in the 
early Universe and for baryogenesis. We hope to return to some of these interesting 
phenomena in the future. 

The example(s) that we analyzed in some detail in this paper were motivated 
by providing a purely four-dimensional description of higher-dimensional fluxes, in 
the spirit of dcconstruction. There arc clearly several possible generalizations of this 
setup. One possibility would be to add flavors to the models discussed in Sect. 3, 
again in the spirit of deconstruction. Note that the net effect of adding such flavors 
is to add new flelds. For example, if we add flavors in the bulk [i.e., at all C/(l) sites 
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except the first and last], the total number of fields for the U{1)^ model would grow 
to 3(n— 1). The naive number of extrema in this case would then be ~ 0(2^"~^) but 
may actually be smaller in practice. There are also new F-term contributions to the 
scalar potential which render some of these vacua stable. 

The field-theoretic framework proposed in this paper also allowed us to explicitly 
calculate the ratio of stable vacua to the total numbers of extrema. We found that 
this ratio is much larger than a naive estimate based on random sign assignments for 
mass eigenvalues would suggest. Wc also found in Sect. 3 that these sorts of results 
— and indeed the entire vacuum structure of the models we considered — are very 
sensitive not only to dimensionful mass parameters, such as FI terms and various 
type of scalar masses (all of which are natural landscape coordinates), but also to 
dimensionless parameters such as the Yukawa and gauge couplings. 

There are clearly more general examples that can arise in string models with D- 
branes. However, these will depend on the particular classes of string models under 
consideration. For example, in heterotic string constructions, there is at most one 
anomalous U{1) factor, with the universal axion mixing a la Stiickelberg [19]. By 
contrast, in a Type I string model (or equivalently a Type II orientifold) , there can 
be a large number of anomalous U{1) factors, each equipped with its own FI term [20]. 

One natural framework to consider consists of taking a large number of f/(l)'s, 
each with its FI term, and all possible charged bi-fundamental fields. If there are n 
different U{1) gauge factors, this would imply n(n+ l)/2 bi-fundamental fields. Such 
a setup would emerge naturally in Type I or Type II string models containing large 
numbers of [/(I) D-branes which either sit at the same point in the compact space or 
intersect each other. In such a geometry, the bi-fundamental fields would emerge as 
open-string excitations stretching from one U{1) brane to another. This setup would 
not only yield a large number of vacua, but would also naturally give rise to various 
new types of superpotential terms. These include cubic, quartic, and higher-order 
interactions, each of which would produce additional F-tcrm contributions to the 
scalar potential. It is easy to imagine that in some region of the parameter space, 
these additional contributions could easily stabilize most of the huge numbers of 
vacua that such models would produce. 

Regardless of the particular field-theoretic model under discussion, it is apparent 
that string models of all sorts generically give rise to multiple U (1) gauge factors. We 
therefore believe that even though the models we presented in this paper are field- 
theoretic in nature, they will emerge naturally in realistic string contexts and should 
be viewed, quite literally, as at least one component of the full string landscape. 
As such, we believe that these types of models should be analyzed further, as they 
undoubtedly play a direct role in affecting the statistical properties of string vacua. 



42 



Acknowledgments 



KRD is supported in part by the US National Science Foundation under 
Grant PHY/0301998, by the US Department of Energy under Grant DE-FG02-04ER- 
41298, and by a Research Innovation Award from Research Corporation. ED is sup- 
ported in part by the CNRS PICS #2530, INTAS grant 03-51-6346, the RTN grants 
MRTN-CT-2004-503369 and MRTN-CT-2004-005104, and by a European Union Ex- 
cellence Grant, MEXT-CT-2003-509661. TG is supported in part by a US Depart- 
ment of Energy grant DE-FG02-94ER-40823, a grant from the Office of the Dean 
of the Graduate School of the University of Minnesota, and a Research Innovation 
Award from Research Corporation. KRD and ED would hke to acknowledge the 
hospitality of the William I. Fine Theoretical Physics Institute at the University of 
Minnesota, and KRD also wishes to acknowledge the hospitality of both CERN and 
the CPhT-Ecolc Polytcchniquc. We also wish to acknowledge the hospitality of the 
Aspen Center for Physics where this work was initiated. 

Appendix 

In this Appendix we present the details behind the n = 3 landscape discussed in 
Sect. 4.1 and sketched in Fig. 11. As discussed in Sect. 4.1, we shall set A = and 
take > 0. The surviving extrema are then in the following classes: 

{0}: These are extrema with all Vi = {(f)i) — 0. The eigenvalues are: 

m^i^i (2), m^iCs (2), (6.1) 

where the number in parentheses corresponds to the multiplicity of the eigenvalues. 
Thus, provided that 

|^i|<m^ and l^sl < , (6.2) 

this vacuum is stable. Clearly, when m = this vacuum is unstable. The D-terms 
are given by 

(^)- {6,0,^3}. (6.3) 
{1}: = ii — m^, V2 = Vs = V4 = 0, with eigenvalues 

0, 2(^1 -m^), 2mM2), ± 6 (2) ■ (6.4) 
The stabihty region is: 

^1 > and |^3| < rh^ . (6.5) 
When rh — this vacuum is no longer stable. The D-terms are given by 

{D)^{m',0,^s} . (6.6) 
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{2}: |v2p = — + 'rn'^),vi = vs = V4 = Q, with eigenvalues 

0, -2(ei + m^), l(Srh'-^^) (2), l(3m^ + - 2^3) (2), m^ + ^a (2). (6.7) 
The stabihty region is: 

Ci < -m^ ^3 > and 2^3 - ^1 - 3m^ < . (6.8) 

When 771 — this vacuum is no longer stable. The D-terms are given by 

(^) = {^(ei-m^),^(6 + m^),?3}. (6.9) 
{3}: \vs\'^ — 1(^3 — m^),vi — V2 — V4 — 0, with eigenvalues 

0, 2(^3 -m^), ^(Sm^ + Ca) (2), ^(377^^ + 2^1-6) (2), - (2). (6.10) 
The stability region is: 

6 < rh^, 6 > rh^, and ^ - 2^ - 3m^ < . (6.11) 
When m = this vacuum is no longer stable. The i^-terms are given by 

{D) = {^uliCs - m%^{^s + rh')} . (6.12) 

{4}: \v4\^ — —^3 — m^, vi — V2 — vs — 0, with eigenvalues 

0, -2(^3 + m'), 2m2 (2), ± 6 (2). (6.13) 
The stabihty region is: 

6 < -m^ and < . (6-14) 
When fh — O this vacuum is no longer stable. The D-terms are given by 

{D) = {^,,0,-m'} . (6.15) 

{12}: There are no solutions in this class. 

{13}: l^ip = — ifi^i 1^3^ = — m^),V2 = = 0, with eigenvalues 

(2), 2(^1 -m^), 2(6 -m'), ^(Sm^-^) (2), ^(Sm' + es) (2). (6.16) 
The stabihty region is 

6 > rn^, and < 6 < 5m^ . (6-17) 
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When m — this vacuum is no longer stable. The D-terms are given by 

(D) ^ {m',\{^s-m%l{^s + rh')} . (6.18) 

{14}: l^ip = ^1 — m^, |v4p = —^3 — m^, V2 — — 0, with eigenvalues 

0(2), 2(6 -m^), -2(^3 + m2), 27h' (A) . (6.19) 
The stability region is 

6 > and ^3 < -m^ . (6.20) 
When m = this vacuum becomes supersymmetric. The D-terms are given by 

(5) = {7fl^ 0, -m^} . (6.21) 

{23}: \V2\' = i(-3m2 - 2^1 + {a), ksP = ^-Sm^ - ^ + 2^3), ^^i = = 0, with 
eigenvalues 

(2), i(6m2±ei±6) (2), 



-4m2 - 26 + 26 ± ^ V3^' + " ^3) + - 6^3 + • (6-22) 
The stability region is 

6<26-3m^ 6>26 + 3m^ and |6 + 6l<6m2- (6-23) 
When m = this vacuum is no longer stable. The D-terms are given by 

0) = {-rh' + ^(6 + 6), ^(6 + 6),m^ + ^(6 + 6)} ■ (6.24) 
{24}: |v2p = |(~Ci ~ "^^)) l'f^4p = "^3 ~ t^i = i'3 = 0, with eigenvalues 
(2), -2(6 + m'), -2(e3 + m2), ^(5m2 + 6) (2), ^(3m2-6) (2). (6.25) 
The stabihty region is 

— 5m^ < 6 < — and 6 < — ■ (6.26) 
When m = this vacuum is no longer stable. The £)-terms are given by 

(D) = {^(6 - rh% ^(6 + rh'), -rh'} . (6.27) 
{34}: There are no solutions in this class. 
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{123}: \vi\'^ = 5 + ^3 - 6m2, l^sp = 6 - Sm^, jvap = ^3 - 37fi2,'f;4 = 0. Some of 
the eigenvalue expressions are quite complicated but one can verify that the stability 
region is: 

+ 6 > and ^3 > 5m^ . (6.28) 

When fh — this vacuum remains stable. The D-terms are given by 

(D) = {m^ 2m^ Sm^} . (6.29) 

{124}: There are no solutions in this class. 
{134}: There are no solutions in this class. 

{234}: |w2p = -Ci-3m^, l^sp = -^i-5m^, 1^4^ = -{i -^3 - 6m^, = 0. Some of 
the eigenvalue expressions are quite complicated but one can verify that the stability 
region is: 

6 + 6 < -6m^ and ^1 < -5m^ . (6.30) 

When fh — this vacuum remains stable. The £)-terms are given by 

{D) = {-3m^ -2m^ -m^} . (6.31) 

{1234}: There are no solutions in this class unless rh — 0. For 7ft = 0, we have a 
line of solutions given by 

= ksl' + Ci, = M-^S: ^2!" = kal' ■ (6.32) 

These solutions are necessarily supersymmetric, with all £)-terms vanishing. 

For < 0, we obtain solutions for the remaining cases: 

{12}: = ^1 — 3m^, \v2\'^ — — 2m^, V3 — V4, — 0, with eigenvalues 

(2), m2+6 (2), 3m'' -^3 (2), -7m^ + 6 ± V25m4 _ q^^^2 + ^2 _ ^5 33) 
This vacuum has no stable regions, and the D-terms are given by 

{D) = {m',2m',^3} . (6.34) 
{34}: |v3p = — 2m^, |v4|^ = —^3 — 2>rh',V\ — V2 — Q.i with eigenvalues 

(2), m'-6 (2), 37ft' + ^1 (2), -7m' - ^3 ± ^aSm^ + 6771^^3 + H ■ (6.35) 

This vacuum has no stable regions, and the D-terms are given by 

{D) = {i,,-2m\-m^} . (6.36) 

{124}: l^ip = —(^3 — m', |w2p = — 3r7i', 1^41' = — 2m',f3 = 0. This vacuum has 
no stable regions, and the D-terms are given by 

{D) = {m^, 2m', -m'} . (6.37) 

{134}: \vi\' = ii — fh', |v3p = —2fh', |u4p = —^3 — 3fh',V2 = 0. This vacuum has 
no stable regions, and the £)-terms are given by 

{D) = {fh', -2fh', -fh'} . (6.38) 
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